Journal of Forestry Research (2009) 20(3): 199-204 
DPI 10.1007/sl 1676-009-0039-6 

RESEARCH PAPER 


Effects of windbreak width in wind direction on wind velocity reduction 


Mulati Yusaiyin, Norio Tanaka 

Graduate School of Science and Engineering, Saitama University, 255 Shimo-okubo, Sakura-ku, Saitama, Saitama, 338-8570, Japan 


Abstract: The variations of drag force acting on the windbreak and the bulk drag coefficients for different windbreak widths were studied 
experimentally in the Eiffel-type non-circulating wind tunnel at the Hydraulic Engineering Laboratory, Saitama University, Japan, to eluci¬ 
date the effects of windbreak width in the wind direction on wind velocity reduction behind a windbreak. The variations of flow field for 
different windbreak widths were studied numerically by using the two-dimensional Reynolds-averaged Navier-Stokes (RANS) equation 
with a k-e turbulence closure model. Results show that the total drag force to wind increased with increasing windbreak width, but the bulk 
drag coefficient decreased slightly. The relationship between the bulk drag coefficient Cd and the windbreak width W and height H can be 
presented by the equation of Cd~kd ( W/H)~ b ( kd, b: constants). The result of the numerical simulation shows that the windbreak width greatly 
affects the location and the value of the minimum wind velocity. The wind velocity decreased by 15%—22% as the windbreak width in¬ 
creased. 
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Introduction 

Windbreaks have been used for centuries to modify wind profiles 
and to defend against wind damage. In the last few decades, 
systematic studies have been considered to determine the effects 
of windbreaks on wind velocity reduction. The commonly used 
descriptor of the internal structure of an artificial windbreak has 
been ‘optical porosity,’ which is defined as the ratio between the 
open area of the windbreak and its total area. Numerous experi¬ 
ments have been conducted with very thin artificial fences or 
screens in laboratory wind tunnels (e.g. Moysey and McPherson 
1966; Raine and Stevenson 1977; Perera 1981) and in field ex¬ 
periments (e.g., Bradley and Mulhearnl983; Wilson 1987). The 
flow field around windbreaks has also been analyzed by numeri¬ 
cal simulations (Hagen et al. 1981; Wilson 1985; Wang and 
Takle 1995). Most recently, Comelis and Gabrieis (2005) meas¬ 
ured the total wind velocity reduction coefficient (1-velocity 
disturbed by the windbreak/velocity in the absence of windbreak) 
of windbreaks with optical porosities (m 2 /m 2 ) ranging from 0.0 
to 1.0, and concluded that a porosity of 0.20-0.35 is optimal in 
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terms of wind velocity reduction considering one, two and three 
rows of windbreak. Thus, the wind velocity reduction effects of 
artificial fences and screens, including the function of windbreak 
porosity, are well understood. 

Under natural conditions, a windbreak forest with a certain 
width will produce aerodynamic effects when wind flows across 
it, distinctly different from that of an artificial wind fence does. 
For thin artificial windbreaks, the porosity is equivalent to the 
optical porosity. For wide natural windbreaks, optical porosity is 
not always the appropriate parameter for expressing the optimal 
condition because optical porosity shows only the 
two-dimensional gap (ratio between open area and total area) 
rather than the three-dimensional spaces (volumetric porosity) 
through which the wind flows across the width of the windbreak. 
Nevertheless, optical porosity has frequently served as a de¬ 
scriptor of natural windbreaks for lack of a practical alternative 
(Heisler and DeWalle 1988). Therefore, most previous studies 
used the drag coefficient for narrow artificial windbreaks (e.g., 
Woodruff et al. 1963; Seginer 1972; Guyot 1978) or measured 
the drag coefficient directly by conducting wind-tunnel experi¬ 
ments on simulated models of narrow and realistic windbreaks 
with different optical porosities (Guan et al. 2003). Only Wang 
and Takle (1996) reported a set of numerical experiments using a 
model system with a fixed porosity windbreak to study the ef¬ 
fects of three dimensionality of the windbreak structure. The 
result showed that the total drag force of the entire windbreak 
varies little with the changes in windbreak width and structure. 
However, windbreak width greatly affects the location of the 
minimum wind velocity, pressure perturbation, and permeability 
of the windbreak. 

The main goal of windbreak design is to find the optimum 
windbreak width that minimizes the loss of land, decreases the 
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cost of planting and maintaining, and yields the optimum protec¬ 
tion of adjacent areas (Wang and Takle 1996). Therefore, the 
objective of this investigation was to study the effects of wind¬ 
break width on wind velocity reduction. The relationship be¬ 
tween the drag coefficient and windbreak width was determined 
by wind tunnel experiments. In addition, the wind flow around 
windbreaks with different widths was also studied using a nu¬ 
merical simulation. 

Materials and methods 

Simulation of surface boundary layer 

Experiments were conducted in the Eiffel-type non-circulating 
wind tunnel at the Hydraulic Engineering Laboratory, Saitama 
University, Japan (Fig. 1). The wind tunnel has a total length of 
11.3 m, and a working section length of 5.0 m with a 
cross-sectional area of 0.5 m width and 0.5 m height. To simulate 
the surface boundary layer, a layer of gravel (diameter = 2-3 cm) 
was uniformly spread on the floor of the first section (length = 
1.0 m) to increase the surface friction and a mean flow velocity 
Uo= 10 m/s was used. The wind profile in the surface boundary 
layer was arranged according to the well-known logarithmic 
function that can be expressed as: 


The friction velocity U t and the surface roughness length z 0 
were estimated from the measurements using equation (1). The 
roughness and friction velocity were 0.002 m and 0.46 m/s, re¬ 
spectively. The boundary layer thickness was 0.10-0.12 m. The 
fluid properties of the simulation should match the turbulent field 
conditions; this can be quantified by the Reynolds number R e as: 



( 2 ) 


where L is the characteristic length (m) (the average of the width 
and height of the wind tunnel working section), u is the mean 
wind velocity, and v is the kinematic viscosity. 

For rectangular pipe, a turbulent regime occurs at Reynolds 
numbers higher than 1 400 (Bagnold 1973). According to this 
value, the flow in our wind tunnel was always fully turbulent at 
velocities greater than 0.049 m/s. Furthermore, the study of 
scaled windbreaks requires that this factor of the wind tunnel 
equals the full-scale factor (Jenson 1958), hence: 
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where H is the windbreak height (m), and subscripts w and / 
denote the wind tunnel and the full-scale conditions, respectively. 
In this study, H/z 0 was set at 500. 


where u z is the wind velocity at height z, U t is the friction ve¬ 
locity, K (=0.4) is the von Karman constant, and z 0 is the sur¬ 
face roughness length. 
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Drag force measurement 


In principle, a windbreak exerts a drag force on the wind field, 
causing net losses of momentum in the incompressible airflow; 
thus, the frontal part of the windbreak reduces the drag force on 
the rear side of the windbreak (sheltering effect) (Raine and Ste¬ 
venson 1977). Although there have been many studies of this 
sheltering effect, most of them investigated thin artificial wind¬ 
breaks and studied the relationship between the optical porosity 
of the windbreak and the drag coefficients. The relationships 
between the drag coefficient and windbreak width have not been 
reported in previous studies. Therefore, we directly measured the 
drag force that changes with windbreak width. A two-axis load 
cell (with stream-wise (X) and transverse (Y) directions, type 
LB-60, SSK Co., Ltd.) that has a resolution of 1/1000 and can 
measure a maximum load of 1.0 Newtons was used to measure 
the drag force on the windbreak model. The drag force per fluid 
mass due to vegetation may be described, as mentioned by Nepf 
(1999), as: 


F = 


force 

mass 


= -C d aU 2 
2 d 


(4) 


Fig. 1 Schematic of vegetation model used in experiments: (a) ex¬ 
perimental apparatus, (b) side view, (c) details of vegetation ar¬ 
rangement (plan view) 


where F is the drag force per fluid mass (N), and C d is the bulk 
drag coefficient, p is the fluid density (kg/m 3 ), U is the equiva¬ 
lent uniform velocity (m/s), and a is the projected plant area unit 
volume (m 2 /m 3 ). If the plants are modeled as cylinders, 
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where N is the number of cylinders unit area (cylinders/nT), 
AS * is the mean space (m) between the cylinders, and d and H 
are the cylinder diameter (m) and height (m), respectively. 

Nepf (1999) defined a dimensionless population density as 
ad= d 2 / AS ~ for a cylindrical model where ad represents 
the fractional volume of the flow domain occupied by plants. 

Windbreak models 

Identical cylinders with diameter (d) of 4.0 mm and height (H) of 
100 mm were used in the experiment in staggered array configu¬ 
rations. The model was installed in the complete cross-section 
width of the wind tunnel from wall to wall perpendicular to the 
wind direction. The experimental apparatus and configuration are 
shown in Table 1 and Fig. 1. 

Table 1. Arrangement of windbreak model and experimental con¬ 
figurations 







Spacing in 

Spacing in 

ad 

Relative windbreak width( W/H) 

cross stream 

direction, S t 

streamwise 

direction, S r 

0.03 

0.20 

0.36 

0.52 

0.68 0.84 

8 d 

4 d 

0.07 

0.14 

0.28 

0.52 

0.88 

4 d 

3 Ad 


Drag forces were measured in every case, and velocity field were measured 
when ad= 0.03, width 0.2 H, 0.84/7. 

Numerical model and approach 

Governing equation 

The simulation was based on the Reynolds-averaged Na- 
vier-Stokes equation (RANS) using the k-e turbulence closure 
model (Launder and Spalding 1974). Four types of equations 
were solved in each case: continuity equation (6), RANS equa¬ 
tion (7), and two turbulence closure equations (8 and 9) for the 
turbulence kinetic energy ( k ) and for the dissipation rate of tur¬ 
bulent kinetic energy (e), respectively. 


du t 

dx, 


= 0 


(6) 
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where U i , U - are the i, j velocity components, respectively, 
X ■ is the j coordinate, p is the air density, G k is the turbulent 


kinetic energy production, and W,W ; - is the Reynolds stress, 
which can be expressed as: 
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(ii) 


where v is the turbulent viscosity, Sy is the Kronecker delta, 
(7 k (= 1.0) and (7 £ (= 1.3) are turbulent Prandt numbers for k 
and e, respectively, and C , C £l , and C £2 are model con¬ 
stants (C = 0.09; C el =1.44; C £l = 1.92; Launder and Spal¬ 
ding 1974). Mean velocity |F| = yj(Uj)~ . 

Boundary conditions 

In this paper, the boundary conditions of the computational do¬ 
main were the inlet boundary conditions used in equation (1). 
The outlet conditions were the homogeneous Neumann condi¬ 
tions u x =v x = 0. At the upper and lower faces, the boundary con¬ 
ditions for k, e and the velocity are wall functions. 

Numerical scheme 

In this model, the Reynolds-averaged Navier-Stokes equation for 
an incompressible fluid is solved in a generalized coordinate 
system. In the present study, the semi-implicit method for pres¬ 
sure-linked equations (SIMPLE) algorithm was used to correct 
the pressure and velocity, and the governing equations were dis¬ 
cretized using the finite volume method. 
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Average velocity 

The averaged wind velocity was calculated as: 


U average =Jj\ H Q U ( Z )dz 


where U(z) is wind velocity at height z. 


( 12 ) 
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Results and discussion 


Validation of numerical simulation 


Relationship between drag force and windbreak width 

The total drag force acting on the windbreak model as a function 
of the relative windbreak width ( W/H) is shown in Fig. 2(a). The 
results show that the total drag force rapidly increased with rela¬ 
tive windbreak widths of less than 0.5 and 0.7 in cases of 
ad= 0.07 and ad= 0.03, respectively. As the relative windbreak 
width increased above those two values, the total drag-force 
increasing rate became lower. The point at where the total 
drag-force increasing rate became less, moved in the windbreak 
downstream direction as the value of ad decreases. This may be 
due to the diminished drag on downstream rows of the wind¬ 
break. This result from the downstream cylinders demonstrated a 
lower impact velocity due to the velocity reduction in the wake 
by upstream cylinders. In addition, the drag force values in¬ 
creased as the windbreak density increased. The bulk drag coef¬ 
ficient decreased slightly as the width increased (Fig. 2b). The 
relationship between the bulk drag coefficient C d and the relative 
windbreak width ( W/H) can be represented by the following 
equation, which is derived from the experimental results shown 
in Fig. 2 as: 

C d =UW/H) h (13) 

where k d , b are constants and depend on the ad value as shown in 
Fig. 2. 



(a) 



Before putting any of these modules into application, key issues 
about them need to be validated and examined to give the best or 
much better performances. To verify the applicability of the tur¬ 
bulence model and the value of the drag coefficient in the present 
study, the wind tunnel experimental data were compared with the 
numerical simulation results. 

Fig. 3 shows comparisons of the simulated and measured ve¬ 
locity vertical profiles located at x = 1 H, 5 H, downstream of the 
windbreak for ad = 0.03, when W/H = 0.2, 0.84. Fig. 4 shows the 
variation of the horizontal velocity profile along the wind direc¬ 
tion at the same height of ziH (= 0.5) for ad= 0.03, when W/H 
=0.2, 0.84, respectively. As shown in these figures, the simulated 
wind velocity results accord well with the experimental values. 



(a) 



(b) 


Fig. 3 Vertical profiles of measured wind velocity in comparison with 
numerical simulation for different streamwise locations: (a) x = 1 H, 
(b) x = 5 H. 
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Fig. 4 Florizontal profile of wind velocity at height z = 0.6 H 

Wind velocity profiles and windbreak width 


Fig. 2 Drag force characteristics for windbreak models with different 
widths: (a) total drag force, (b) bulk drag coefficient 
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The velocity field was calculated by changing the relative wind¬ 
break widths from 0.2 to 8 using the numerical simulation model, 
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to study the effect of windbreak width on the wind velocity re¬ 
duction behind the windbreak. Characteristics of wind velocity 
reduction behind the windbreak can be expressed by many pa¬ 
rameters. The most commonly used are shelter length (Z 2 o), over 
which wind velocity in the downstream is reduced by 20%, and 
the wind velocity minimum value ( U min ) and its location (X,„ /n ) as 
shown in Fig. 5. 



Fig. S Definition of characteristic parameters that explain the shelter 
effect 



(a) 



(b) 

Fig. 6 Horizontal profile of wind velocity for different values of 
windbreak width at z = (1.5//: (a) for «</=0.03, (b) for «rf=0.07 

The effects of changes in windbreak width on the wind veloc¬ 
ity horizontal profiles are shown in Fig. 6 for different windbreak 
widths at height 0.5 H , and in Fig. 7 for the same widths at 
height 0.15 H , when ad= 0.03, 0.07, respectively. As the wind¬ 
break width increased, the location of the minimum wind veloc¬ 
ity moved toward the windbreak and may move into the 

wider windbreak (negative X mi „). With increasing windbreak 
width, the wind velocity first decreased until it reached its mini¬ 
mum value, and then increased again. When ad= 0.03, the rela¬ 
tive windbreak width value of 7 had the lowest wind velocity 
value, while for ad= 0.07, the relative windbreak width value of 5 
had the lowest wind velocity value. These results are similar to 


those of Wang and Takle (1996). However, at the height of 0.15 
H, the downstream velocity recovered quickly as the windbreak 
width increased, and the minimum wind velocity (U mln ) was 
higher than that at the height of 0.5 H. This result differs from 
that of Wang and Takle (1996). 

Fig. 8 shows the relationships between L 20 and windbreak 
width for ad= 0.03 and 0.07, respectively. The results in Fig. 8 
show that L 20 significantly increased until the relative windbreak 
width was equal to 1 for ad~0.01 or to 3 for ad= 0.03. Further¬ 
more, it decreased slightly as the width increased. It was found 
that the windbreak widths greatly affected L 20 values because 
windbreaks with higher density and narrower width have larger 
L 2 o values. 



(a) 



(b) 


Fig.7 Wind velocity horizontal profiles for different values of wind¬ 
break width z = 0.15//: (a) for ad=0.03, (b) for ad=0.07 



Windbreak width (H) 


Fig. 8 Values of L w for different windbreak widths. 

The average relative velocity values with different windbreak 
widths in experiments at streamwise locations x = 10 H and 15 H 
are shown in Fig. 9 (a) and (b), respectively. The results show 
that at location x = 10 H the average velocity clearly decreased 
until the relative windbreak width was equal to 3 for ad= 0.03, 
and to 1 for ad= 0.07, and then slightly increased as the width 
increased. The windbreak with the relative width of 1 for 
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ad=0.01 had the lowest average wind velocity at location x = 10 
H. Fig. 9(b) clearly shows that the average velocity decreased 
when the relative windbreak width was less than or equal to 1, 
and slightly increased as the width increased for both ad values 
(ad= 0.03 and 0.07). 



(a) 



Windbreak width 
(b) 

Fig. 9 Relative average velocities for different windbreak widths: (a) 
location x = 10 //, (b) location x = 15 H 

Summary and conclusions 

The variations of the drag force and the bulk drag coefficient 
values with different windbreak widths were experimentally 
determined to investigate the effects of windbreak width on the 
wind velocity profile behind a windbreak. The variations in the 
flow with different windbreak widths were also studied numeri¬ 
cally. The relationship between the bulk drag force coefficient C d 
and the relative windbreak width (W/H) can be presented by the 
function of Cj=k d ( W/H) b . The drag force increased with in¬ 
creased windbreak width, but the increase was not linear because 
of the sheltering effect. The windbreak width greatly affected the 
minimum wind velocity and location behind and downstream 
from the windbreak. 
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